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Ground-state properties of the spin-1 two- leg antiferromagnetic ladder are investigated precisely 
by means of the quantum Monte Carlo method. It is found that the correlation length along the 
chains and the spin gap both remain finite regardless of the strength of interchain coupling, i.e., 
the Haldane state and the spin-1 dimer state are connected smoothly without any quantum phase 
transitions between them. We propose a plaquette-singlet solid state, which qualitatively describes 
the ground state of the spin-1 ladder quite well, and also a corresponding topological hidden order 
parameter. It is shown numerically that the new hidden order parameter remains finite up to the 
dimer limit, though the conventional string order defined on each chain vanishes immediately when 
infinitesimal interchain coupling is introduced. 
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I. INTRODUCTION 

Quantum spin-ladder systems have been studied theo- 
retically and experimentally over the last decade as ma- 
terials with a novel spin-gap state, as well as by their rel- 
evance to the high-temperature superconductivity.EI Es- 
pecially, the two-leg ladder Heisenberg antiferromagnet, 
which is defined by the Hamiltonian: 

H — J {Si^i • Si^i+i -I- 82,1 • S2,i+i} 

i 

+/^^Si,, -82.^, (1) 

i 

has been studied most extensively. Here, SaA is the spin 
operator at site i on the a-th chain (a = 1,2), and the 
intrachain and interchain coupling constants are denoted 
by J and K , respectively. In the following, we restrict our 
attention only to the case in which the intrachain cou- 
pling is antiferromagnetic (J > 0). On the other hand, 
the interchain coupling constant, K, can be either posi- 
tive (antiferromagnetic) or negative (ferromagnetic). 

At = 0, the system consists of two decoupled anti- 
ferromagnetic Heisenberg chains. In this case, it is well 
known that the ground-state properties can be classified 
into two universality classes depending on the parity of 
2S. Here S is the spin size. In the case where S is 
a half-odd integer, the ground state is critical, i.e., the 
system has gapless low-lying excitation and the antifer- 
romagnetic correlation function along the chain decays in 
an algebraic way as the distance increases. On the other 
hand, it is conjectured by Haldanea that the antiferro- 
magnetic Heisenberg chain of integer spins has a finite 
excitation gap above its unique ground state, and the 
correlation function decays exponentially with a finite 
correlation length. Its ground-state properties can be 
understood quite well frorruthe viewpoint of the valence- 
bond solid (VBS) picturejj in which the ground state 



is essentially represented as direct products of spin-i 
dimers (AKLT state). In addition, a topological order 
parameter characterizing the AKLT state, as well as the 
Haldane state, ^so-called the string order parameter, has 
been proposed.Q The validity of Haldane's conjecture has 
been confirmed preiasely for S — 1, 2, and 3 by several 
numerical methods Hla 

Introduction of non-zero interchain coupling, K, is 
known to drastically change the ground state, at least for 
the spin-i case.ElFor small K, either antiferromagnetic or 
ferromagnetic, it immediately opens_a.spin gap of Od^l) 
with some logarithmic corrections.EJ That is, K — 
is the special point at which there occurs a quantum 
second-order phase transition between the dimer phase 
{K > 0) and the spin-1 Haldane phase (K < 0). Again, 
from the viewpoint of VBS picture, one can understand 
this phase transition as a global rearrangement of dimer 
pattern. For larger half-odd-integer spins (5 = |, |, • • •), 
the criticality a,t K = should be essentially the same as 
in the spin-^ case. 

In the case of integer-spin chains, on the other hand, ef- 
fects of interchain coupling have been known little so far. 
Recently, Senechal and Allen studied tiic spin-1 ladder by 
mapping it to the nonlinear a modeltll and also by the 
bosonization technique.E3 They found that in contrast 
to the spin-i case, small interchain coupling reduces the 
magnitude of the spin gap in both of antiferromagnetic 
and ferromagnetic regimes. In addition, their analyses as 
well as their complemental Monte Carlo calculation sug- 
gest that there is no critical point between the Haldane 
and the spin-1 dimer phases. This may seem paradoxical 
since these two phases have apparently different dimer 
patterns from each other. 

In this paper, we present the results of our exten- 
sive quantum Monte Carlo simulation on the spin-1 lad- 
der. After reviewing details of our simulation using the 
efficient continuous- time loop algorithm in Sec. we 
present our numerical data on the uniform susceptibil- 
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ity, staggered susceptibility, antiferromagnetic correla- 
tion length, etc. in Sec. Ill, which convincingly demon- 
strate the continuity of the two limiting case {K = 
and K = oo), and thus support, [the conjecture by the 

Ml 



(a) 



In Sec. IV 



we pro- 



previous analytical approaches. 
pose a plaquette- singlet solid state, which is constructed 
as products of local singlet states of four S — ^ spins. 
The Haldane state and the spin-1 dimer state are natu- 
rally included as special limits. In addition, we propose 
a kind of hidden order parameter, which can detect the 
topological hidden order exsisting in the plaquette-singlet 
solid state. We show numerically that the new hidden or- 
der parameter we propose remains finite in the whole pa- 
rameter range, < X < oo, while the conventional strin, 
order parameter vanishes except at = 0. In Sec 
we consider in turn the case where the interchain cou- 
pling is ferromagnetic {K < 0), and show that the spin-1 
Haldane state and that of S" = 2 also continue to each 
other without any singularity on the way to the other. 
We give a summary of our results and some discussions 
in the final section. 



II. QUANTUM MONTE CARLO METHOD 
A. Continuous-time loop algorithm for spin-1 system 

The-, secently-developed continuous- time loop algo- 
rithmll3il3 is one of the most efficient methods for simu- 
lating quantum spin systems. It is a variant of the world- 
line Monte Carlo method, which is based on the path- 
integral representation by means of the Suzuki- Trotter 
discretization.ll3 However, the continuous-time loop alm^ 
rithm works directly in the imaginary-time continuum,t£l 
and thus is completely free from the systematic error in 
the Suzuki- Trotter discretization. In addition, the cor- 
relation between successive spin configurations is greatly 
reduced, sometimes by orders of magnutide, since it flips 
effectively clusters of spins, or loops, whose linear sizes 
correspond directly to the length scale of relevant spin 
fluctuations. The algorithm has alreadyr-been applied to 
various spin systems with great success .1111 

The Hamiltonian we consider is given by Eq. (|l]) with 
S = 1. The linear size along the chain is denoted by L, 
and we adopt periodic boundary conditions in this direc- 
tion, i.e., Sa^i+L = ^a,i for a = 1 and 2. In order to apply 
the continuous-time loop algorithm to the present spin- 
1 system, first we represent the spin-1 Hamiltonian in 
terms of subspins. In this representation, each spin-1 op- 
erator in Ea. (|l|) is decomposed into a sum of two spin-i 
operators .Ea Simultaneously, each bond of strength J (or 
K) is transformed into four bonds of the same strength 
connecting subspins. The lattices before and after the 
subspin transformation are shown in Fig. |l|. Note that in 
order to recover dimensions of the original spin-1 Hilbert 
space (3^^), one needs to introduce a set of projection 
operators, each of which acts on a pair of subspins and 
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FIG. 1. (a) Original spin-1 ladder, and (b) equivalent sys- 
tem represented by 5 = i spins (subspins). The spin-1 lad- 
der of length L {2L spins and 31/ bonds) is mapped onto the 
spin-i system on a lattice of 4L subspins and 12L bonds. 
Each oval in (b) denotes a pair of subspins which is sym- 
metrized by being applied special boundary conditions in the 
imaginary-time direction. 



projects out the state with 5' = (Fig. After trans- 
formed into a path-integral representation, the projection 
operators are converted to speciaLboundary conditions 
in the imaginary-time directionp'EEl for each pair of sub- 
spins, the total is required to be conserved across the 
imaginary-time boundary. 

For the mamied., system, the spin-i continuous-time 
loop algorithmE^nlj can be applied without any modifi- 
cation except that we need to introduce additional graphs 
and labeling rules for the boundaries in the imaginary- 
time direction .EIIIj We use the multi-cluster variant of the 
loop algorithm. The resulting algorithm is found to work 
quite well as the same as the original algorithm developed 
for S* = -i; the integrated auto-correlation time for the 



physical quantities we measure remains of order unity, 
and there is observed no significant sign of its growth even 
in the largest system in the present simulation {L = 1024 
and T/(J + K) = 2/1024 ~ 0.00195). 



B. Physical quantities 

The physical quantities of interest can be measured by 
using the corresponding subspin representations. First, 
for later convenience, we introduce the imaginary-time 
dynamical correlation function, C^{x, r), and its Fourier 
transform, C^{k,uj). The former is defined explicitly in 
the path-integral representation by 

x{Sl,^+,{t + T)±S2,^+At + T)}) , (2) 



where (3 is the inverse temperature (= l/T). The spin 
configuration at site i on the a-th chain and at imaginary 
time T (0 < T < /3) is denoted by Sa^iir), which takes -1, 
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0, or 1. The bracket (• • •) in Eq. (|^) denote the average 
over Monte Carlo steps (MCS). In the present subspin 
representation, Sa,i{T) is simply given by a sum of 5*^ of 
two subspins at (a,i,T). 

In terms of the imaginary-time dynamic structure fac- 
tor, C^{k,uj), the uniform susceptibility and the stag- 
gered susceptibility are simply given by 



and 



X = 2L/3C+(0,0) 



Xs = 2L/3(7-(^,0), 



(3) 



(4) 



respectively. We also calculate the static structure factor 
at momentum tt as 



^(E(-1)'^"''('5m(0)- 52,(0)) 



S{n) 

x(5i,,(0)-52,,(0))^ 

In order to calculate the correlation leng; 
chains, we use the second-moment methodi 




27r/i,0) 



(5) 

along the 
(6) 



Similarly, the spin gap, which is defined as the inverse of 
the correlation length in the imaginary-time direction, is 
measured by 



A- 



C^(7r,0) 
C'±(7r,2V/3) 



(7) 



In Eqs. and (|^), we take the minus (plus) sign for 
if > (K < 0). Although the above second-moment 
estimates suffer from systematic error due to the exis- 
tence of subdominant decaying modes in the correlation 
function, it should be sufficiently small (at if = the 
systematic error for the spin gap is known to be about 
0.29?tl), and thus we expect that it would be irrelevant to 
the following discussions. We will also present our results 
for the string order parameteiu and a new hidden order 
parameter in Sec. IV. Their explicit definitions will be 
given later. 

In practice, all the physical quantities we will show in 
the following, including the hidden order parameters, are 
measured by using so-called improved estimators. For 
example, the staggered susceptibility for if > is simply 
represented as the sum of squared length of each loop, 
divided by 8Lf3. Similarly, the imaginary-time dynamic 
structure factor can be measured directly as 



(8) 



where the integration is performed along a closed path 
on each loop, and the summation runs over all loops. 



TABLE I. Convergence of physical quantities at J = 0.7 
and K = 0.3. The temperature is taken as T/{J + K) = 2/ L 
for each L. The figure in parentheses denotes the statistical 
error (2(j) in the last digit. There are observed no significant 
difi'erences between the data with L = 256, 512, and 1024. 



L 
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(J + if) Xs 




A/(J + A) 
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3 X lO'' 


16.48(3) 


4.95(1) 


0.3863(6) 


4.57(1) 


16 


6 X 10^ 


49.72(8) 


7.74(1) 


0.2120(3) 


8.50(1) 


32 


7 X 10^ 


138.0(2) 


11.39(2) 


0.1218(2) 


14.86(2) 


64 


1 X 10'' 


315.9(4) 


14.97(2) 


0.0793(1) 


22.81(3) 


128 


1 X lO'' 


468.7(8) 


16.23(2) 


0.06602(9) 


27.37(4) 


256 


7 X 10^ 


489.8(5) 


16.19(2) 


0.06487(7) 


27.91(4) 


512 


3 X 10^ 


490.3(3) 


16.17(2) 


0.06477(8) 


27.93(3) 


1024 


2 X lO ^' 


490.1(2) 


16.17(2) 


0.06476(4) 


27.92(2) 



C. Taking the thermodynamic limit and the 
zero-temperature limit 

In the present method, the system size. L, and the 
temperature, T, are restricted to be finitcE^ Since we are 
mainly interested in the ground-state properties of the 
infinite lattice, a proper extrapolation scheme for taking 
the thermodynamic limit (L — > oo) as well as the zero- 
temperature one (T — > 0) is required. In the present 
study, we adopt the following strategy, which is the same 
as was used in Ref. ^ for the estimation of the Haldane 
gap of the antiferromagnetic Heisenberg chain with S = 
1, 2, and 3. 

For each parameter set {J,K), we start with a small 
lattice at relatively high temperature, e.g., L — 8 and 
T/{J + K) = 0.25. Then, we increase the system size ex- 
ponentially step by step, like L = 16, 32, 64, 128, • • •. Si- 
multaneously, the temperature is decreased so as to keep 
LT/{J -I- K) = LI{J -I- if)/3 constant. In other words, 
the aspect ratio of the (l-l-l)-dimensional space is kept 
unchanged for all i's. If the system is gapless, the finite- 
size scaling holds; the correlation length in the real-space 
direction, S^{L,T), and that in the imaginary-time direc- 
tion^ A"i(L,r), both would grow being proportional to 
L.E3 Simultaneously, the other physical quantities, such 
as the staggered susceptibility, should exhibit power-law 
behavior with some exponents depending on their own 
anomalous dimensions. 

On the other hand, if the system is gapful (this is the 
case for the present system as we will see below), there 
exist finite intrinsic correlation lengths, ^ and A~^. As 
long as the system size and the inverse temperature are 
smaller enough than these intrinsic correlation lengths, 
the critical behavior above mentioned is still observed. 
However, once both of L and /3 exceed them enough, 
£,{L,T) and A(L,T), as well as other physical quanti- 
ties, no longer exhibit system-size dependence. Strictly 
speaking, the systematic error due to the finiteness of 
the system decreases exponentially as the system size in- 
creases, and it becomes much smaller than the statistical 
error due to the finiteness of MCS. 
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In Table |, we show the system-size (and tempera- 
ture) dependence of the physical quantities at (J, K) = 
(0.7,0.3). As seen clearly, the data with L > 256 (and 
with T/{J+K) < 2/256 = 0.0078125) exhibit no system- 
size dependence. Thus, in this case, one can safely con- 
clude that the system is gapful and also that the phys- 
ical quantities obtained are those of the infinite lattice 
at zero temperature besides the statistical error. Empir- 
ically, we find that L/^{L, T) > 6 with /3A(L, T) > 6 is a 
reasonable condition to guarantee the convergence in the 
present numerical accuracy. 

One of the advantages of the present scheme is that it 
depends on no numerical extrapolation techniques, such 
as least-squares fitting, the Shanks transform, etc. Final 
results are simply obtained from those of the largest sys- 
tem at the lowest temperature in the simulation. There- 
fore, this method is quite stable and the error estimation 
is also quite reliable. It should be emphasized that the 
final results are not affected at all by the value chosen for 
the aspect ratio, L/{J + K)(3 {—2 in the present case). 
However, if one chooses a too small or too large value, 
the physical memory of the computer system might be 
exhausted before reaching the thermodynamic limit or 
the zero-temperature one. 

In what follows, we will mainly present the data with 
L = 256 and r/(J K) = 2/256 unless otherwise 
noted. Measurement of physical quantities is performed 
for 5 X 10^ - 1 X IQ^ MCS after 10^ MCS for ther- 
malization. Typically, simulation of this system requires 
about 7MB of physical memory, and 1 MCS takes about 
0.33 sec. on a single CPU of SGI 2800 (MIPS R12000 
400MHz). This system size is somewhat exaggerated 
for certain sets of parameters, {J,K). However, in such 
cases, the system can be considered as Af statistically- 
independent samples simulated in parallel, where J\f ~ 
(L/^) X (A/3). Therefore, we gain better statistics pro- 
portional to V^, which completely compensates-for the 
growth of CPU time (~ LP for large L and /3).Ej Thus, 
we loose nothing besides the memory requirement. This 
is already manifested in the figures presented in Table |. 



III. NUMERICAL RESULTS 

A. Parameterization 

The ground state of the present system is parametrized 
by the ratio of the interchain coupling constant to the 
intrachain one, x = K/J. Hereafter, we mainly use i?, 
which is defined by 



0.3 



R = 



K 



J +1X1 



1 



(9) 



Since we consider only the antiferromagnetic intrachain 
coupling (J > 0), -1 < 7? < 1. At i? = 0, the 
system consists of two independent antiferromagnetic 
chains (spin-1 Haldane chains). On the other hand, at 
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FIG. 2. Temperature dependence of the uniform suscep- 
tibility for 7? = (squares), 0.2 (circles), 0.3 (upward tri- 
angles), 0.4 (diamonds), 0.6 (downward triangles), and 0.8 
(pentagons). The exact result for R — 1 (Eq. (|lO|)) is plotted 
by a solid line. In the inset, the same data are plotted against 
the logarithm of T. 



i? = 1, it is decoupled into dimers sitting on each rung. 
The limit R ~* ~1 corresponds to a single spin-2 an- 
tiferromagnetic chain. We also introduce the reduced 
temperature, T ~ T/{J + K), reduced susceptibilities, 
X = {J + K)x and Xs = {J + K)xs, and the reduced gap, 
A = A/(J -1- K). In this section and Sec. IV, we con- 
sider only the case where < i? < 1, i.e., the interchain 
coupling is antiferromagnetic. 



B. Uniform susceptibility 

Before investigating the zero-temperature properties 
of the spin-1 ladder, first we briefly discuss finite- 
temperature behavior of the uniform susceptibility, which 
gives a rough profile on the excitation spectrum of the 
system. In Fig. ||, x is plotted as a function of T for 
R ^ 0, 0.2, 0.3, 0.4, 0.6, 0.8, and 1. At i? = 1, i.e., 
J = 0, since the system consists of independent dimers, 
the exact form of the uniform susceptibility is easily ob- 
tained as 

, _ 1 exp(-l/f ) + 5exp(-3/T') 
^~ f l + 3exp(-l/f) + 5exp(-3/f) 



-1/f) forf<l. 



(10) 



At i? = 0, X is also known to vanish at low temper aiures 
exponentially as exp(-A/f ) with A ~ 0.41050(2)0 be- 
sides a prefactor of some powers of T. As seen clearly in 
Fig. H, the temperature dependence of the uniform sus- 
ceptibility does not depend on the value of R strongly; it 
has a very broad peak around T ~ 1, and decreases quite 
rapidly at lower temperatures. It suggests spin-singlet 
ground state regardless of the value of R. 

However, it should be noted that at temperatures lower 
than T ~ 0.4, the uniform susceptibility is not a mono- 
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FIG. 3. 7?-dependence of the staggered susceptibihty (open 
squares) and the static structure factor at momentum vr (solid 
circles). The statistical error of each data point is much 
smaller than the symbol sizes. 
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FIG. 4. i?-dependence of the inverse correlation length 
along the chains (solid circles) and the spin gap (open 
squares). The statistical error of each data point is much 
smaller than the symbol sizes. 



tonic function of R. It is greatly enhanced by orders of 
magnitude around i? = 0.3, in comparison with those at 
i? = and 1. It indicates that in the intermediate region 
of i?, the spin gap is strongly suppressed, and on the 
other hand long-range antifcrromagnetic fluctuations are 
enhanced, due to the competition between the intrachain 
and the interchain antiferromagnetic couplings. 

In addition, in the temperature profile of the uniform 
susceptibility for 0.2 < R < 0.6, one can see a clear 
s/iowWer structure at T ~ 0.1 (see also the inset ofFig.||), 
which indicates existence of some additional anomalies in 
the excitation spectrum of the system. 



C. Quantities at zero temperature 

The temperature dependences of the staggered sus- 
ceptibility, Xs, and the static structure factor, ^(Tr), are 
found to be qualitatively very similar to those of the sin- 
gle Haldane chain {R = 0); they grow rapidly around 
T ~ 1 and are saturated to finite values at low tem- 
peratures, though the saturation temperature strongly 
depends on R (see the i?-dependence of A shown below). 
As for S{t:), in addition, it makes a weak peak before 
saturated to the zero-temperature value, which indicates 
a short-range antiferromagnetic order. 

In Fig. H, the zero-temperature values of Xs and S{tt) 
are plotted as a function of R. The value of Xs at i? = is 
consistent with that obtained in the previous work, Xs — 
18.4048(7)0. On the other hand, at i? = 1, they coincide 
the exact values, Xs = f and S'(7r) = |, respectively, 
within the statistical error. One sees in Fig. ^ that they 
are smooth functions of i?, and thus there is no indication 
of singularities in the whole range of R. 

The non-existence of phase transitions is also con- 
firmed by the i?-dependence of the spin gap, A, and the 
inverse correlation length, (Fig. §). These results 
convincingly support the conjecture made in the previous 



analytic studies. It should be recall ed th at by using 
the method we explained in detail in Sec. II C , the conver- 



gence to the thermodynamic and the zero-temperature 
limits of the data has been checked for all the value of R 
we simulated. Therefore, the data can be identified with 
those at L = CX3 and T = besides the statistical errors, 
which are much smaller than the symbol sizes in Figs. |3| 
and ||. 

Although there are no singularities between R = and 
1, long-range antiferromagnetic fiuctuations are greatly 
enhanced in the intermediate region of R. This is consis- 
tent with the temperature de pende nce of the uniform sus- 
ceptibility presented in Sec. [II B. Especially, note that 



all the physical quantities we calculated have its maxi- 
mum (or minimum) at i? ~ 0.3 (see also Table |). They 
could be compared with those of the spin-2 antiferromag- 
netic Heisenberg chain, i.e., Xs = 1164.0(2), ^ = 49.49(1), 
and A = 0.08917(4)i 



IV. PLAQUETTE-SINGLET SOLID STATE AND 
HIDDEN ORDER PARAMETER 

A. Breakdown of AKLT picture 

As we have already mentioned in Sec. ||, the ground 
state of the spin-1 antiferromagnetic chaiQ.is understood 
quite well by means of the VBS picture. El Actually, the 
AKLT state, which is the exact ground state of the so- 
called AKLT modeljj shares many common properties, 
such as spin-rotation symmetry, finite correlation length, 
etc., with the ground state of the spin-1 Heisenberg chain. 
These two states are believed to belong to the same uni- 
versality class with each other. 

The AKLT state is essentially constructed as di- 
rect products of spin-i dimers sitting on each bond 
(Fig. 11(a)). On each site, two edge S = ^ spins are 
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FIG. 5. Schematic picture of plaquette-singlet solid state, 
(a) AKLT state at 7? = 0. (b) Spin-1 dimer state at 7? = 1 
(c) Plaquette-singlet solid state for < i? < 1. 



symmetrized to form an 5 = 1 spin. An important fea- 
ture of this state is that a sum of in any interval 
can take only 1, 0, or -1. This fact is immediately fol- 
lowed by the existence of a topological hidden order; if 
one removes spins at 5^ = state, the sequence of of 
the remaining spins has a perfect antiferromagnetic-like 
ordering (1,-1, 1,-1, •••): although the original one has 
no 'true' antiferromagnetic long-range order at all. 

The above topological order in the AKLT state can 
be detected, quantitatively by means of the string order 
parameterp which is defined by 



{O2 



hm {02{t,j)) 

-j\-*oo 



in terms of the string correlation operator: 



02(*,j) = -^fexp 



(11) 



(12) 



whose expected value is | = 0.444 • • • for any \i — j\ > 
1. The string order parameter ( [Tl] ) is also finite in the 
spin-1 Heisenberg chain. Its value is estimated to be 
0.3743(1) in the thermodynamic limit (i? = in Fig. H). 
In practice, the string order parameter of a finite chain 
of length L is calculated as 



2 L 



\Y.^02{^,^ + L/2)). 



(13) 



Although the string correlation function (|T^) is a non- 
local quantity, one can construct an improved estimator 
even for it. For details, see Ref. 

In the presence of the interchain coupling, however, 
the above AKLT picture breaks down immediately. As 
shown in Fig. |6[ the string order parameter, (02)l, de- 
fined on one of the two chains, decreases quite rapidly as 
R increases. Furthermore, in contrast to the other quan- 
tities shown before, the value of {02)l exhibits strong 
system-size dependence for i? > 0. Actually, as shown in 



the inset of Fig. ^, one finds that these finite-size data 
scale exponentially quite well as 



{O 



21L 



fiR log"L) 



(14) 



with a = 2.5, where f{x) is a scaling function, and it 
vanishes for x — > 00. This strongly suggests that the 
string order parameter ( pi] ) is essentially singular at i? = 
0, and is vanished by infinitesimal interchain coupling, 
though we have no account for the value of the exponent 
a at the moment. 



B. Plaquette-singlet solid state 

As we have seen above, the AKLT picture (Fig. || (a)) 
breaks down immediately for R > 0. Actually, in the 
R = 1 limit, the ground state is represented schematically 
by a pattern of spin-^ dimers, in which two dimers are 
sitting on each rung (Fig. I (b)). There is no overlap 
between these two VBS states. 

However, if we focus our attention on a plaquette con- 
taining four spins at {a,x) = (l,i), (l,z -I- 1), (2,i), and 
(2, i + 1), it is found that these two states share a remark- 
able feature, i.e., the sum of of these four spins can 
take only 0, ±1, or ±2, and its absolute value never ex- 
ceeds 2 in both cases. In other words, four S = ^ spins, 
each of which belongs to an 5 = 1 spin at one of four cor- 
ners of a plaquette, form an 5 = state, though its fine 
structure is quite different from each other. This obser- 
vation leads us to a proposal of the following generalized 
state, which connects these two VBS states smoothly. 

First, we define a local state of plaquette consisting of 
four S — ^ spins, which is expressed explicitly as 

|V(0)), = ag[cose{\ T)i,,| Di.rn - I T)i..+i} 

X {|T)2,,U)2,»+1 - U)2,,:|T)2,»+l} 

+ sin0{|t)i,.U)2,.-U)i,.|T)2.} (15) 
X {|T>i.«+iU>2,,:+l - I 1)2,1+1} 

where ag is the overall normalization factor (a^"^ = 
\/4 -f 2 sin 29). This state is constructed as a superpo- 
sition of two S — states. The first term of r.h.s. in 
Eq. ( p^ is the product of the dimers sitting on legs, and 
the second one does that of the rung dimers. The pa- 
rameter 9 {0 < 9 < controls the proportion of these 
two constituents. Note that the wave function (|l5| ) is 
the exact ground state of the following four-body spin-i 
Hamiltonian: 

"Hp — -^{(Ti,i ■ Cri,i+1 + 0-2,4 • 0'2,j+l) 



■ 02 



(16) 



where Jp > 0, > 0, and cr^,! denotes the spin-iPauli 
operator. In this case, the parameter 9 in Eq. (15) can 
be expressed explicitly as a function of Xp = Kpjjp as 
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FIG. 6. _R-dependence of the hidden order parameters, 
(O2) (open symbols) and ^ (O4) (solid circles). The statisti- 
cal error of each data point is much smaller than the symbol 
sizes. The horizontal dashed line denotes the exact value (|) 



of ^ {O4) in the dimer limit (7? — 1). A scaling plot for (O2) 
with exponent a = 2.5 is also shown in the inset. 



d = arctan 



(17) 



and the gap above the ground state never closes through- 
out while one varies Xp from (6* = 0) to cx) (0 = -1), as 
we will see in Sec. VL 

By using the local singlet state (|T5|), we construct a 
wave function of the spin-1 ladder: 



(18) 



where Va^i is a projection operator acting on two S = =^ 
spins at site (a, i). The schematic picture of this state is 
presented in Fig. || (c). We refer to it as the plaquette- 
singlet solid state. One sees that by varying the parame- 
ter 9 from to -I, it connects smoothly the AKLT state 
at i? = (Fig. ^ (a)) and the spin-1 dimer state at i? = 1 
(Fig. H (b)). We expect the ground state of the present 
system can be described well by the plaquette-singlet 
solid state with 9 tuned for each specific value of i?.E3 

Furthermore, it is possible to define a four-body string 
correlation operator. 



•E 

k=i+l 



'SI 



S: 



2,k) 



QZ QZ 
'^l,j'^2,j J 



(19) 



which characterizes the plaquette-singlet solid state (18). 
It is easy to prove that the expectation value of 04(1, j) 
remains finite for any value of 9 for the plaquette-singlet 
solid state (H). Especially, (OiiiJ)) = | at S = f (i? = 
1). It is interesting to see that ©4(1, j) can be expressed 
as a product of the conventional string correlations (|lT 
defined on each chain, 02{i, j, a): 



(20) 
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FIG. 7. Demonstration of topological hidden order in the 
plaquette-singlet solid state, (a) Spin configuration on each 
plaquette. Arrows denote the configuration (±i) of S = ^ 
spins. On each plaquette, the sum of configurations of four 
spins is restricted to be zero, (b) Spin configuration of cor- 
responding spin-1 ladder, (c) Sequences {Si} and {%} cal- 
culated from the configuration (b) as Si — Sf^i + S^^i and 
Ti = iSl j, respectively, (d) Sequence of Ti's obtained by 
being removed elements satisfying both of 7i = and \Si\ 7^ 1, 
and replaced 1 and -1 by (1,0) and (-1,0), respectively, in 
which all positive (negative) Ti's sit on one (another) sublat- 
tice of the sequence. 



Therefore, in the decoupled-chain case (i? = 0) we have 

{04{l,j))^{02{l,j))\ (21) 

where we omit the index a, since the expectation value 
of the string correlation operator 02(1, j; a) does not de- 
pend on a. For R > 0, such a simple relation does not 
hold, and in general they take different values with each 
other. 

In Fig. ^, we show the i?-dependence of our new hidden 
order parameter. 



(a 



4 L 



i + L/2)) 



(22) 



calculated for the present model. We plot {OilL in- 
stead of {04)l itself to demonstrate the relation ( pl| ) at 
R = Q. In contrast to the conventional string order pa- 
rameter, {02)l, which vanishes immediately for R > 0, 
the new string order parameter, (04)^, is found to be a 
smooth function of R, and remains finite up to the dimer 
limit (R = 1). 

It should be emphasized that although the long-range 
spin fluctuations aregreatly enhanced round R = 0.3 
as discussed in Sec. HI, (04)l remains remarkably large 
even in this region. This strongly supports that the 
ground state of the spin-1 ladder is actually described 
quite well by the plaquette-singlet solid state (^8|) in the 
whole range of R. 
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R 

FIG. 8. i?-dependence of the spin gap (open squares) and 
the inverse correlation length along the chain (solid circles) in 
the ferromagnetic {K < 0) case. The dashed lines denote the 
respective values in the R—>—l limit (Ref. ||). 



0.41 - 

0.4 - slope = 3.75(3) / \slope = -3.72(3) - 

0.39 - 
0.38 - 
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0.37 ^ 1 ' 1 

-0.01 -0.005 0.005 0.01 

K/J 

FIG. 9. if/ J-dependence of the spin gap, A/ J, in the anti- 
ferromagnetic {K > 0) and ferromagnetic {K < 0) cases. The 
dotted straight lines are obtained by least-squares fitting for 
< K/J < 0.025 and -0.025 < K/J < 0, respectively. 



Before closing this section, it might be worth not- 
ing that the plaquette-singlet solid state (Eq. ( p^ ) and 
Fig. ^ (c)) has the following topological antiferromagnetic 
long-range order, though its definition is slightly compli- 
cated than its counterpart in the AKLT state for the 
single spin-1 chain; first, for a given spin configuration 
{Saj} {ct — 1,2 and i = 1, 2, 3, • • •), define two sequences, 
(5i,52,53, • • •) and (Ti,7^,7^, • • •), where 5i's and 7^'s 
are calculated as Si — ^+S2 , and % = SI ji"! ^, respec- 
tively. Next, choose a nearest pair of non-vanishing T^'s, 
say Tj and 7^. If the number of S^s satisfying = 1 in 
the interval j < ^ < fc is even (odd) , then Tj and 7^ has 
a same (opposite) sign. In other words, in the sequence 
of Ti's, if one removes elements satisfying both of 7^ = 
and \Si\ ^ 1, and replaces 1 and -1 by (1,0) and (-1,0), 
respectively, then one finds the resulting sequence {%} 
has an antiferromagnetic long-range order regarding non- 
zero elements, i.e., all positive (negative) T^'s sit on one 
(another) sublattice of the sequence. A demonstration of 
this topological order is presented in Fig. 0. 

V. FERROMAGNETIC INTERCHAIN 
COUPLING 

So far, we consider only the case where the interchain 
coupling is antiferromagnetic (i? > 0). In this section, we 
consider in turn the ferromagnetic case. In the i? — > — 1 
limit, two S ~ 1 spins connected by an infinitely-strong 
ferromagnetic rung bond form S = 2. Thus, two intra- 
chain antiferromagnetic coupling of strength J in Eq. (Q) 
are transformed into one bond of strength J = J/2 con- 
necting two S* = 2 spins. The spin-2 antiferromagnetic 
chain has a finite spin gap. A/ J — 0.08917(4), which is 
much smaller than that of the spin-1 |ehain, and also a 
longer correlation length, ^ = 49.49(1) S 

The spin-2 Haldane state is also explained by a VBS 
picture, in which two dimers sitting on each bond. There- 



fore, one would expect naturally that the two indepen- 
dent spin-1 Haldane state at i? = and the spin-2 one 
(i? — > —1) is connected smoothly without any singulari- 
ties between them. Indeed, as shown in Fig. ^ the spin 
gap. A/ J, as well as the correlation length, decreases 
monotonically as R does, and seems to converge smoothly 
to the values of the spin-2 Haldane chain. Note that the 
spin gap presented in Fig. ^ is not the reduced gap, A. 
Since A vanishes as ~ l/|Ar| for R 0, we normalize 
the gap as A/ J, which remains finite in the i? — > — 1 
limit. At i? = -0.8 and -0.9, we need L = 512 and 
T = 0.0009766 to obtain the zero-temperature values in 
the thermodynamic limit. 

Although the decoupled-chain point {K = 0) is not a 
critical point, the spin gap is still non- analytic at this 
point. In Fig. ||, we show the X/J-dependence of A/ J 
near K = {\K/J\ < 0.01) both in the antiferromagnetic 
and ferromagnetic regimes. One sees a clear cusp just at 
K = 0. In both regimes. A/ J decreases linearly as \K/J\ 
increases. Furthermore, we find the absolute value of its 
slope is the same within the error on both sides (3.72(3) 
ioT K > and 3.75(3) for K < 0). We should remark 
this is completely consistentp-with the conjecture by the 
previous bosonization study.t3 

VI. SUMMARY AND DISCUSSIONS 

In this paper, we investigated precisely the ground- 
state properties of the spin-1 antiferromagnetic Heisen- 
berg two-leg ladder by means of the quantum Monte 
Carlo simulation. We found that the system is gapful 
regardless of the strength of interchain coupling, that is, 
the Haldane state in the decoupled chains and the spin- 
1 dimer state are connected smoothly and there is no 
quantum phase transition between them. We conclude 
that the behavior of the spin gap as a function of R we 
observed, including the cusp at i? = 0, is completely 
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consistent with the recent analytic, studies based on the 
mapping t|e-jthe nonUnear a modellij and the bosonization 
technique Jl3 

Although there is no quantum phase transition in the 
present system, the spin gap is greatly suppressed, and at 
the same time the long-range antiferromagnetic fluctua- 
tions are enhanced in the intermediate region (i? ~ 0.3), 
e.g., A/( J + K) = 0.06476(4) and (J K)x^ = 490.1(2) 
at i? = 0.3. This explains the reason why the spin-1 
ladder with i? ~ 0.3 exhibits an antiferromagnetic long- 
range order, if one introduces quite small rfc; 0.004 J) 
antiferromagnetic coupling between laddersBa It should 
be noted that the critical inter-ladder coupling is smaller 
by an order of magnitude than the spin gap. This situ- 
ation is the same as the critical interchain coupling, J^, 
of the two-diinensional array of spin-1 chains, J'^j J — 
0.043648(8),y while A/J = 0.41050(2).! The precise 
ground-state phase diagram of these two-dimensional 
spin-1 Heisenberg antiferromagnets has been reported in 
detail in Ref. |2^. 

The AKLT picture, which works quite well for the sin- 
gle spin-1 chain, was found to break down immediately 
in the presence of interchain coupling. We proposed the 
plaquette-singlet solid state to describe the ground state 
qualitatively for the whole range of R. The plaquette- 
singlet solid state has a topological hidden order. The 
hidden order parameter (©4), which characterizes the 
plaquette-singlet solid state, was shown to be finite up 
to the dimer limit (i? = 1). This strongly supports that 
the ground state of the spin-1 ladder is described quite 
well by the present plaquette-singlet solid picture. 

It is quite interesting to see that the local Hamiltonian 
of spin-i plaquette, Eq. (|l^), itself can reproduce quali- 
tatively correct behavior of the spin gap observed in the 
present spin-1 ladder. The Hamiltonian (16) has a sin- 
glet ground state regardless of the value of Xp = A'p/ Jp 
(—00 < Xp < 00). The Xp-dependence of the spin gap 
can be written explicitly as follows: 



A, 



for Xp > 



1 - Xd 



Xp for Xp < 0. 



(23) 



One can find easily the following features: (i) There al- 
ways exists a finite gap for — 00 < Xp < cxd. (ii) For 
Xp > (ATp > 0), A/( Jp -I- ATp) has a minimum at a finite 
value of Xp. (iii) On the other hand, for Xp < (ATp < 0), 
A/Jp is a monotonically-decreasing function of jxpl, and 
converges to a finite value at Xp = —00. (iv) Finally, at 
Xp = 0, A/Jp has a cusp, and it decreases linearly with 
the same slope on both sides. The last behavior is simply 
due to the crossing of the second- and third-lowest eigen- 
values. All of these features are qualitatively the same as 
the present spin-1 ladder. 

The plaquette-singlet solid picture proposed in the 
present paper might be generalized to wider classes of the 
spin-1 Heisenberg ladder. Indeed, for the ferromagnetic 
interchain coupling case discussed in Sec. ^ our hidden 



order parameter (O4) was found to remain at a finite 
value (~ 0.089) in the limit K —00. Furthermore, the 
present plaquette-singlet solid picture might be applied 
even in the presence of bond alternation (forced dimeriza- 
tion) in the intrachain coupling, i.e., Ji = (1 + (— 1)'(S) J, 
where S denotes the strength of dimerization. We have 
observed that there is no singularity for the whole range 
of (5 (0 < (5 < 1), and (O4) remains finite up to the 
decoupled-plaquette limit {6 = 1). Note that as for the 
decoupled dimerized chain {K — 0), the (5-dependence of 
(O4) is qualitatively different from that icf^-K > 0. There 
exists a critical point at Sc = 0.26001 (4) |23 and the con- 
ventional string order parameter (O2) vanishes for 6 > 6c. 
Since the relation ( |2l] ) holds for AT = 0, the new hidden 
order parameter (O4) also vanishes for S > Sc. The ex- 
istence of non-zero hidden order parameter (O4) in the 
cases with ferromagnetic couplings and with alternating 
bonds implies their ground states are also well described 
by products of local singlet states similar to the present 
plaquette-singlet solid state, though the structure of lo- 
cal singlet state should be quite different from the present 
one. 

Very recently a new nitroxide material, abbreviated 
as BIP-TENO, has been synthesized, and itS|-mag- 
netic properties have been investigated precisely.E§Ej The 
tetraradical BIP-TENO molecule consists of two pairs 
of ferromagnetically-coupled S — ^ spins, and relatively 
weak antiferromagnetic coupling exists between the pairs. 
The crystalline state of BIP-TENO is thus expected to 
be described effectively by a spin-1 antiferromagnetic lad- 
der of present interest. From the high-temperature be- 
havior of the uniform susceptibility, Katoh et at esti- 
mated the strength of effective couplings and jConcluded 
that J ~ 50K and AT ~ 42K, i.e., R ~ 0.46.EI A clear 
spin-gapped behavior of the uniform susceptibility ob- 
served in the experiment is consistent with the present 
results. The excitation gap A is ajlsp estimated to be 
15. 6K from the magnetization curve.E^ Unfortunately, its 
magnitude is more than twice as large as the present re- 
sult (A/(J + K) = 0.07392(6) at i? = 0.45. See Fig. |). 

Interestingly, the shoulder in the temperature depen- 
dence of the uniform susceptibilitjLXFig- §) has also been 
observed even in the real material.ES This anomaly in the 
excitation spectrum might be related with the ^-plateau 
observed in the magnetization process of BIP-TENO. o 
However, it is beyond the scope of the present study, and 
remains as a future problem. 
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